Effective spin-wave action for ordered Heisenberg antiferromagnets in a magnetic field 
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We derive the effective long-wavelength Euclidean action for the antiferromagnetic spin-waves of 
ordered quantum antiferromagnets subject to a uniform magnetic field. We point out that the mag- 
netic field dependence of the spin-wave dispersion predicted by the usual 0(3)-quantum nonlinear 
sigma model disagrees with spin-wave theory. We argue that the nonlinear sigma model does not 
take into account all relevant spin-wave interactions and derive a modified effective action for the 
long-wavelength spin-waves which contains an additional quartic interaction. At zero temperature 
the corresponding vertex is relevant in the renormalization group sense below three dimensions. 
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I. INTRODUCTION 

Recently, two of us derived a new effective field theory 1 
for transverse spin fluctuations in ordered quantum an- 
tiferromagnets. Our approach is based on the Holstein- 
Primakoff transformation, 2 which maps the original spin 
Hamiltonian onto a bosonic many-body Hamiltonian. 
However, rather than working directly with the canoni- 
cal Holstcin-Primakoff bosons, we express them in terms 
of Hermitian field operators representing staggered and 
uniform transverse spin fluctuations^ The main advan- 
tages of this procedure are: (a) the suppression of the 
effective interaction between antiferromagnetic magnons 
at long wavelengths is manifest; (b) the relation between 
Holstein-Primakoff bosons and the fields of the nonlinear 
sigma model^ (NLSM) representing the transverse stag- 
gered spin fluctuation is made precise; (c) the zero wave 
vector modes of the uniform and staggered magnetization 
can be easily treated within this approach, 3 which is im- 
portant in the study of finite magnets; (d) the transverse 
ferromagnetic fluctuations may be eliminated to yield an 
effective theory involving only transverse antiferromag- 
netic fluctuations. This is most easily done within a path 
integral approach; the resulting effective action coincides 
with the NLSM at the Gaussian level. Although the in- 
teraction terms of the effective action differ from those 
of the NLSM, the renormalization group flows of both 
theories agree to one-loop order i 

Here, we extend the Hermitian operator approach to 
antiferromagnets in a uniform magnetic field. In this 
case, the differences between this approach and the 
NLSM become more substantial. Already at the Gaus- 
sian level the two theories differ: while the Hermitian 
operator approach reproduces by construction the cor- 
rect spin-wave dispersion known from spin-wave theory, 
the NLSM does not. More severely, we show below that 
a quartic interaction term present in the Hermitian op- 
erator approach is absent in the NLSM. 

We shall constrain the analysis here to a Heisenberg 
model on a hypercubic lattice with nearest neighbor an- 
tiferromagnetic coupling J > 0. It is however straight- 
forward to extend the results to general bipartite lat- 



tices. This would also cover a recently synthesized metal- 
organic 5 = 5/2 quantum antiferromagnet on a distorted 
honeycomb lattice^ which is bipartite. In presence of a 
uniform external field h (measured in units of energy), 
the Hamiltonian is given by 



H = jY,Si-Sj-^2hSi 



(1) 



Here, Si are spin operators normalized such that Sf = 
S(S+ 1), and the first sum is over nearest neighbor pairs 
(ij) of a hypercubic lattice with lattice constant a. 

Our main interest is a description of the low-energy and 
long- wavelength properties of the Heisenberg model l[T]l. 
A standard approach to this problem is the 0(3) quan- 
tum NLSM whose effective action describes only the rele- 
vant staggered spin fluctuations. According to Refs. 00 , 
a magnetic field can be taken into account within the 
NLSM approach by replacing the derivative d T with re- 
spect to the imaginary time r with a covariant derivative 
(or Lie derivative, see Ref. 



d T — > d T — ih x 



(2) 



With this minimal coupling the NLSM in a uniform mag- 
netic field has the form 



'NLSM 



[«] = £ 



r(d T Sl-ihxft) 2 , (3) 



where the unit vector 0(t, r) represents the slowly fluc- 
tuating staggered magnetization, po and cq are the spin 
stiffness and spin-wave velocity at temperature T = 
1/(3 — 0, and 9 M = d/dr^ is the spatial derivative in 
direction p = 1, ...,D. Numerical values of po and 
Co must be computed microscopically, e. g. by means 
of a l/S'-expansioni* To leading order in 1/S, one has 
pa fa JS 2 a 2 ~ D and c fa 2D 1 / 2 JSa. For later reference, 
we note that the classical uniform transverse susceptibil- 
ity is xo = Po/ c o with xo ~ (ADJa )^ 1 to leading order 
in 1/S. The classical uniform magnetization per volume 
is M = xoh- 
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FIG. 1: Spin configuration (Si) = Srhi in the classical ground 
state of a two-sublattice antiferromagnet subject to a uniform 
magnetic field h = he x , where Ti belongs to sublattice A 
Here mo = h/(4:DJS) and no = 



II. SPIN- WAVE DISPERSION 



In absence of a magnetic field, the NLSM predicts 
a doubly degenerate antiferromagnetic spin-wave mode 
with long- wavelength dispersion Ep. — Crj|fe|, in agree- 
ment with linear spin-wave theory,* Rather surpris- 
ingly, however, in presence of a magnetic field the long- 
wavelength dispersion derived from the NLSM in Eq. © 
disagrees with spin-wave theory. To see this, let us re- 
call the spin-wave dispersion predicted by linear spin- 
wave theory. In this approach the spin operators are ex- 
panded in deviations from the classical ground state spin 
configuration, using either the Holstein-PrimakofS or 
the Dyson-Malee^ transformation. The classical ground 
state of a quantum antiferromagnet in a uniform external 
field is canted, where the staggered magnetization is per- 
pendicular to the magnetic field, as shown in Fig. To 
leading order in l/S, spin- wave theor y^' 1 ?? 11 predicts two 
transverse spin- wave modes, one gapless and one gapped, 
with long- wavelength dispersions (see Eqs. (|17I20|I below) 



E k+ w [h 2 + elk 2 } 1 / 2 , E k _ 
where the spin- wave velocities are 



c+ = c \ 1 - 3mg , c_ = c a \ 1 - mg 



(4) 



(5) 



Here, mo = |too| is the length of the normalized classical 
uniform magnetization 



m 



D M 



h 



S ADJS ' 



(6) 



which is related to the classical canting angle i?o be- 
tween the direction of the staggered magnetization and 



the local magnetic moments by mo = sin^o- I n general, 
the classical canting angle is renormalized by quantum 
fluctuations £1H 

Let us compare Eqs. Ill o f) with the dispersion rela- 
tion obtained from the NLSM. Denoting by e a the unit 
vectors in the fixed directions a — x, y, z and assuming 
h = he x to point along the x-axis, we write 



n = me x + Uj_e y + [i - nf - njj 



2U/2, 



(7) 



and expand the action Snlsm[0] to quadratic order in 
the fluctuations II m and n_L transverse to the direction 
of the staggered magnetization. One easily finds from 
Eq. © that the spin- wave mode ELi polarized parallel to 
the magnetic field is gapped with energy dispersion 



E H = [h 2 + c 2 k 2 Y' 2 



(8) 



while the spin- wave mode II j_ polarized perpendicular to 
the magnetic field has a gapless linear dispersion, 



E k ± = co\k\ 



(9) 



The spin-wave spectrum of the NLSM given in Eqs. (JSJ 
and iJHJ does not reproduce the correct magnetic field- 
dependence of the spin-wave spectrum in Eqs. illl- r >l) ob- 
tained via conventional spin-wave theory. Obviously, the 
above minimal coupling J2J does not account for order 
(h/J) 2 corrections to the dispersion, even though the 
qualitative aspects of the dispersion are correctly de- 
scribed. 



III. EFFECTIVE ACTION 

In Ref. we have developed a new method to de- 
rive the effective action for staggered transverse spin- 
fluctuations. We here sketch the main steps of the deriva- 
tion adapted to a quantum antiferromagnet subject to 
a uniform magnetic field^ Starting point is the rep- 
resentation of spin operators in terms of canonical bo- 
son operators hi and b\ using the Holstein-Primakoff 
transformation 2 * Let us denote by rhj = (Si)/S the 
directions of the spins in the true ground state and intro- 
duce a local right-handed orthogonal triad of unit vec- 



tors e 



(i) 



(2) 



As explained in Ref. 01 there is a 
U(l) gauge freedom in the choice of the transverse vec- 
Defining the spherical basis vectors 

ef + ipef , p = ±, we express the components of the 
spin operator Si in terms of canonical boson operators bi 
and b J as follows™ 



tors e,- 1 ^ and 



,( 2 ) 



2 h 



with 



si = s-m 

S+ = {2S) 1/2 [1 



rii = b\bi , 
~2s) b 



(10) 

(11a) 
(lib) 



3 



and 5/~ = (S^~) . The Heisenberg model £Q) can then 
be written as a bosonic many-body Hamiltonian^ We 
expand around the classical ground state configuration 
(see Fig. [TJ 



and choose 



J 1 ) 



.(2) 



m e z 



(12) 



(13) 



where Q = 1 ( d 
n = y/l - m 2 = 
part of this Hamiltonian, wc first perform a Fourier trans- 
formation in the sublattice basis 



= -1) for n € A (n € B). Here, 
cosi^o- To diagonalize the quadratic 



_ J (2/N) 1 / 2 J2 k e lk ^A k , neA 
(2/N)^ 2 j: k e lk ^B k , neB 



(14) 



where the momentum sums are over the reduced (an- 
tiferromagnetic) Brillouin zone. After introducing the 
symmetric and antisymmetric combinations C ka = 
2~ 1 l 2 [A k + aB k ], a = ±, we apply the Bogoliubov trans- 
formation 



Cfe 



— ko 



Uka 
-OVkG 



where 



Uka = 



Vka = 



Uka 



T-O^k + Cka 



rl 



2e ka 
ml"f k - tka 



^ka 

^ k 

— ko 



1/2 
1/2 



with 



('kc 



= [(1 + aml 1 kf - («o7fc) 



21I/2 



(15) 

(16a) 
(16b) 

(17) 



Here, j k = D 1 cos(fc ■ a M ), where a M , \i — 1 . . . D, 
arc the D primitive lattice vectors of the hypercubic lat- 
tice. The quadratic part of the effective boson Hamilto- 
nian assumes then the form of non-interacting harmonic 
oscillators 



H2 — ^ Eka 



k.a 



(18) 



where the operators ^> k(T satisfy the canonical bosonic 
commutation relations, 



(19) 



and the energy dispersion of the two spin- wave modes are 
given by 



Eka = 2D J Seka ■ 



(20) 



In the limit of long wavelengths Eq. 1|2(J[1 reduces to the 
energy dispersions given in Eqs. 1 Hot . 



To derive the long-wavelength effective action for the 
staggered transverse spin fluctuations, we need the pre- 
cise relation between the Bogoliubov quasi-particle oper- 
ators ^ka and the field operators Uka representing the 
transverse fluctuations of the staggered magnetization. 
In a Euclidean path integral approach, these operators 
correspond to continuum fields Uka (t) which are analogs 
of the fields n fe ||(r) and n fc± (r) of the NLSM. The rela- 
tion between these two different parameterizations of the 
spin fluctuations is simply 1 

*fe+ = -i{xo/2VE k+ ) 1/2 [E k+ nk++iXo 1 ®k+] ,(21a) 
*fc- = (xo/2KE fc _) 1/2 [Ek-tlk- +^x ( 7 1 $fc-] , (21b) 

where V = Na D is the volume of the system. One easily 
verifies that the pairs II CT , <3E><x satisfy commutation rela- 
tions of canonically conjugate bosonic field operators, 



[n fc(T , $k>a>} — iV5k.-k'5a,a> 



(22) 



The field operators Uka and <5> ka have a simple interpre- 
tation in terms of transverse staggered and uniform spin 
components. The Fourier transformed spin operators on 
the A and B sublattice are given by 



Sa/bm = (2/iV) 1/2 e~ ik ' ri Si 



(23) 



riGA/B 



Using Eqs. i|10ll3fl we find that, to leading order in 1/5 
and mo, the transverse staggered components of the spins 
are given byi* 

S$ k = ^(4V S$ k ) « (S/a D )N-^Tlk- , (24a) 
sSk = ^iSf k + Sf M )^-{S/a D )N-^n k+ .{2Ab) 



In contrast, the transverse components of the uniform 
magnetization are related to the operators <& kcr , 



,(1) _ J_^ Q (i) 
/2 

o(2) _ 1 /o(2) _ o(2) x 

~~ y/2 °B,k) 



S k 1} = ^(S^ k + S^ k )^N-y^ k+ , (25a) 

N- 1/2 § k - • (25b) 



In terms of these operators, our quadratic spin-wave 
Hamiltonian l|18l) can be written as 



H 2 = 



2V ^ [ Xo 1 



^-ka^ka+XnE 2 ka ti-katika ■ (26) 



The effective Euclidean action 5 G ff[n cr ] for the stag- 
gered spin fluctuations can now be obtained by writing 
the partition function as a phase space path integral over 
the fields H ka (r) and $/ co .(t) associated with the above 
operators and subsequently integrating over the "Iv-fields 
which represent gapped ferromagnetic fluctuations. At 



4 



the level of a Gaussian approximation, we obtain in this 
way S c ft « S-g with 

s$[n CT ] = Jg- ^> 2 + £L)n-^n KCT . (27) 

Here, we have combined momenta k and bosonic Mat- 
subara frequencies u> n in a composite label K = (fe, iw n ). 
We have further defined 

n Ka = / dr e^« T n feCT (r) . (28) 
Jo 

At long wavelengths the effective action (|27fl has the same 
form as the corresponding Gaussian part of the action of 
the NLSM. However, in contrast to the NLSM Q, our 
effective action (|27|l has the correct spin- wave dispersion, 
even for short wavelengths. 

We have calculated the leading corrections to the effec- 
tive action S^sp^] for staggered spin fluctuations aris- 
ing from spin- wave interactions in the Holstein-Primakoff 
approach^ Keeping only interaction terms which be- 
come relevant below three dimensions and thus omitting 
terms which are marginal in D = 1, wc find in the con- 
tinuum limit 

s off [n CT ] = -^^( w 2+ c 2fc 2 + r CT )n_^n A v 

- ixoh dr J d D rTP + d T Ti._ 

- d D rn 2 + (H 2 + +n 2 _). (29) 

Here, the c a are the Holstein-Primakoff results for the 
spin-wave velocities, which for large S are given in 
Eq. (JSJ. To leading orders in 1/5, the values of the 
gap parameters are r + = h 2 and r_ = 0. For symme- 
try reasons, the true spin- wave spectrum of the n_-mode 



must remain gapless at vanishing wave- vector (Goldstone 
mode), 15 while the spin- wave gap of the n + -mode is not 
renormalizecU 5 *^. In a renormalization group analysis, 
this requires fine tuning such that at the fixed point these 
conditions are met. While the cubic term in Eq. (|29|l cor- 
responds precisely to the cubic (Berry phase) term in the 
NLSM © , the quartic interaction between the field com- 
ponents is absent if the magnetic field is included in the 
NLSM by the minimal coupling 

IV. SUMMARY AND CONCLUSION 

In this work we have used the Holstein-Primakoff 
transformation to derive an effective action for staggered 
transverse spin fluctuations of quantum Heisenberg an- 
tiferromagnets in uniform magnetic fields. Our effective 
action contains an additional quartic interaction between 
the field components which is not contained in the NLSM. 
It is easy to see that at zero temperature the quartic 
interaction vertex in Eq. I|29|l is relevant in the renor- 
malization group sense for dimensions D < 3, so that 
we expect that it generates singularities in perturbation 
theory^ We conclude that the NLSM given in Eq. 
does not contain all relevant interactions in the ordered 
phase below three dimensions. The reason for this is 
that not all approximations that are usually made in the 
derivation of the NLSM 6 arc justified in the presence of a 
uniform magnetic field. In particular, it is not justified to 
neglect Umklapp scattering processes between transverse 
and longitudinal spin fluctuations involving momentum 
transfers across the boundary of the magnetic Brillouin 
zone. As an alternative to Eq. Ij29(l it should be possible 
to analyse the long-wavelength staggered spin fluctua- 
tions within a phenomenological Ginzburg-Landau model 
with a quartic interaction termii since this model treats 
the longitudinal modes as independent degrees of free- 
dom. 
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